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Three topics in the Schwinger model* 
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1. We compare Monte Carlo results with analytic predictions for the fermion condensate, in the massive one- 
flavour Schwinger model. 2. We illustrate on the Schwinger model how to facilitate transitions between topological 
sectors by a simple reweighting method. 3. We discuss exact, non-perturbative improvement of the gauge sector. 

1. FERMION CONDENSATE 

In the massless (one-flavour) Schwinger model, 
the fermion condensate is 



= e^/27r = 0.283466, 



(1) 



in units of n = g/ ^/tt. For non-zero fermion mass 
m, no exact expression is known, but there are 
some results from analytical approximations and 
lattice simulations [|l|. No clear agreement has 
been reached, though. 

Here we report on a simulation using the R- 
algorithm |^ with staggered fermions, in which 
extrapolations to zero step-size, to zero lattice 
spacing and to infinite volume are performed. We 
focus on two values of m/ fi: 0.33 and 1.00. 

Fig. shows the approach to the continuum 
limit for a given physical volume: the condensate 
is plotted against 1//3 for a series of simulations 
on {L/a) X {2L/a) lattices {L/a = 6, 8,..., 16), 
keeping am\/]3n = 0.33 and {L/a)/y/j3 — 4^2 
fixed. The condensate for the free theory is sub- 
tracted (in two different ways, see caption) in or- 
der to cancel the UV divergence present in the 
massive case. A linear extrapolation to the con- 
tinuum limit, accounting for 0{l/f3) ^ C(a^) 
perturbative scaling violations, works well and 
is seen to be very important, given the large ef- 
fect. After additional extrapolations to infinite 
volume and zero MD step-size, we find for the 




'Contribution to Lattice '97 by A.J. van der Sijs 
(arjanOscsc.ethz.ch). ETH-SCSC preprint TR-97-12. 



0.2 0.4 



Figure 1. Fermion condensate against 1/(3, 
with extrapolation, for fixed physical volume and 
fermion mass. The two sets of points correspond 
to subtraction of the free condensate in the same 
finite volume (solid circles) or in infinite volume 
(open circles). The four squares correspond to 
2x2 times as large volumes. 



(subtracted) condensate: 

(-i^^)sub/M = 0.141(5) (m//i-0.33), (2) 

and similarly for the larger mass: 

(-V5V')sub^ = 0.084(3) (m/fi^im). (3) 

In fig. 1^ these data are compared with a re- 
cent analytical result by HosotaniQ expected 

•^We thank Y. Hosotani for discussions on this topic. 
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Figure 2. The fermion condensate in infinite vol- 
ume, in the continuum limit (squares) compared 
with the analytical small-m result of Ref. |^ (solid 
curve). The two asterisks correspond to our re- 
analysis of data from Ref. fl. 



to hold in the small-m region. For comparison, 
we have also plotted our reanalysis {i.e., subtrac- 
tion of the free condensate and subsequent linear 
extrapolation to 1/(3 — 0) of data from Ref. Q 
for the non-compact model on a 64 x 64 lattice. 

2. SAMPLING TOPOLOGICAL 
SECTORS AT HIGH /3 

A correct sampling of the different topologi- 
cal sectors is essential for a correct determina- 
tion of quantities like the fermionic condensate 
in the Schwinger model. In the usual impor- 
tance sampling algorithms, however, transitions 
between sectors are suppressed by a factor of or- 
der exp(2/3) in the continuum limit, implying ex- 
ponential critical slowing down for the topological 
autocorrelation time. 

We propose the following reweighting technique 
to overcome this problem. We add to the action 
the quantity 



AS* 



E 
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(4) 



and compensate for this by including a factor 
exp(/3AS') in the observables. Here xq is the site 



Figure 3. Topological charge history (in units of 
M.D. time) for a simulation on an 18 x 36 lattice 
at /? = 10.125, using reweighting. For compar- 
ison, without reweighting we observed only one 
transition in 200 000 units of M.D. time. 



at which the plaquette angle is closest to tt 
of all the plaquettes in the configuration under 
consideration. This procedure increases the topo- 
logical transition rate by strongly enhancing the 
probability that the plaquette at xq goes through 
TT. Note that Xq may vary from one configura- 
tion to the next, but the essential point is that 
the prescription assigns a unique action to each 
configuration. It can be shown that this leads to 
reversible molecular dynamics, ensuring detailed 
balance. The computational overhead is negligi- 
ble. 

The example in Fig. ^ shows that this method 
drastically improves sampling of the topological 
sectors (ergodicity). 

The parameters a and can be used for op- 
timization. For example, they can be adjusted 
such that the distribution (histogram) of 0^1/ (a;o) 
is roughly flat over most of the interval [— 7r,7r]. 
On large lattices, 9^v{xo) gets closer to tt, such 
that smaller a and 0o suffice for the same transi- 
tion rate, and the bias due to Eq. (^ gets smaller. 

This reweighting procedure enhances a local 
lattice artefact, which acts like a "saddle point" 
between topological sectors in the lattice formu- 
lation; in the continuum theory, these sectors 
are completely disconnected. This observation 
may guide a simple generalization to non-abelian 
gauge theories in 4 dimensions. 
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3. NON-PERTURBATIVE 
IMPROVEMENT 

Two-dimensional U(l) lattice gauge theory 
with the standard plaquette action is exactly solv- 
able. In particular, in the infinite-volume limit 
the expectation value of any rectangular Wilson 
loop of charge q satisfies an exact area law with 
string tension 
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a'^(Tq{f3) = In 
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where Io.q{l3) are the Bessel functions. For the 
non-compact action, on the other hand (cf. also 
the Manton action), one finds (Tq — q'^/2f3 ("per- 
turbative scaling") exactly. We use this observa- 
tion as a guideline for non-perturbative improve- 
ment^ of the compact plaquette action: tune the 
coefficient of the improvement term such that the 
string tension equals /2(3. 

We consider the addition of an adjoint (squared 
plaquette) term to the standard action: 



A5 = ^c(/3) 



1 — cos ( 



1 — cos 20p 
12 



(7) 



In this case, the string tension is given by Eq. (H) 
with /„ replaced with A„, defined as 



A„(/3/,/3a) 



d9 



cos(n6l)e''^™"^+'^''™"2e^ (g) 



with = ^ (1 c(/3)/3), /3a = -/3c(/3)/12. The 
\n are calculated numerically. 

In Fig. ^ we compare various choices for c(/3). 
By taking c = 1, the 0^ term in the action, 
and hence the term in c^Oq (|6|), are can- 
celled (perturbative improvement): the slope of 
the relative "error" in the string tension (fig. 
^) increases from 1 to 2. The error is seen 
to be further reduced when this coefficient is 
tadpole-improved, by taking c(/3) = l/(n(/3)) 
self-consistently. Non-perturbative improvement 
amounts to tuning c(/3) such that a^Cg = <z^/2/3 

^This is different from tlie usual improvement where the 
expansion of the plaquette angle in terms of the continuum 
gauge field plays an essential role. 
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Figure 4. c(/3) for the improvement schemes 
described in the text (a) and relative deviation of 
a^CTi from the value 1/2/3 (6). For the dash-dotted 
curve in (&), the horizontal axis is l/2/3eff. 



for all 13 (for a given q) . Note that a very good re- 
sult is obtained by simply using the unimproved 
(c = 0) Wilson action but expressing the results 
in terms of the effective coupling /3cff = /3 (□(/3)) 
suggested by the tadpole scheme. 
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